Angular dependence of coercivity in magnetic nanotubes 
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The nucleation field for infinite magnetic nanotubes, in the case of a magnetic field applied parallel 
to the long axis of the tubes, is calculated as a function of their geometric parameters and compared 
with those produced inside the pores of anodic alumina membranes by atomic layer deposition. We 
also extended this result to the case of an angular dependence. We observed a transition from curling- 
mode rotation to coherent-mode rotation as a function of the angle in which the external magnetic 
field is applied. Finally, we observed that the internal radii of the tubes favors the magnetization 
curling reversal. 

PACS numbers: 75.75.-|-a, 75.10.-b 



I. INTRODUCTION 

Since the discovery of carbon nanotubes by lijima in 
1991, [i| intense attention has been paid to hollow tubu- 
lar nanostructures because of their particular significance 
for prospective applications. More recently, magnetic 
nanotubes have been grown 0, 0, U, |^, 0] motivating in- 
tense research in the field. Technological applications of 
such systems require a deep knowledge and characteriza- 
tion of their magnetic behavior. For example, changes in 
the internal radii are expected to strongly affect the mag- 
netization reversal mechanism [3] and thereby the over- 
all magnetic response. 0, [1] The nature of the magnetic 
tubes may be suitable for applications in biotechnology, 
where magnetic nanostructures with low density, which 
can float in solutions, are very desirable. 3 

Coercivity is one of the most important properties of 
magnetic materials for many present and future appli- 
cations of permanent magnets/magnetic materials, mag- 
netic recording, and spin electronics and, therefore, the 
understanding of magnetization reversal mechanisms is a 
permanent challenge for researchers involved in studying 
the properties of these materials. Recently, Landeros et 
al. [7| found that, when a magnetic field is applied paral- 
lel to the axis of a tube, the curling reversal mode is the 
dominant magnetization reversal mechanism for tubes 
with radii greater than 30 nm. However, the angular 
dependence of the coercivity in magnetic nanotubes has 
not been studied yet, in spite of many works on this topic 
comprising nanowires. [s, [l3| In this work we calculate 
the angular dependence of the coercivity of Ni nanotubes, 
assuming that the reversal of the magnetization occurs 
by means of one of two possible modes: magnetization- 
curling mode [V) and coherent-rotation mode (C). 

Geometrically, tubes are characterized by their exter- 
nal and internal radii, R and a, respectively, and length 
L. It is convenient to define the ratio (3 = a/R, so that 
f3 = represents a solid cylinder and (3 close to 1 corre- 
sponds to a very narrow tube. Besides, we consider an 
external magnetic field Ha applied in a direction defined 



by 6*0, with do the angle of the applied field with the tube 
axis, as illustrate in Fig. 1. In this paper we present an 
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FIG. 1: Geometrical parameters of a magnetic nanotube. 

analytical model about the switching modes and fields of 
infinite extended magnetic nanotubes in dependence of 
the orientation of the magnetic field versus the nanotube 
axis. Additionally, experimental data for the switching 
field of high-aspect ratio Ni nanotubes, when the mag- 
netic field is applied parallel to the tube axis, will be 
compared with this micromagnetic model. 

II. EXPERIMENTAL DETAILS AND RESULTS 

The high-aspect ratio Ni nanotubes were produced in 
porous membranes with pore diameters of 180, 220 and 
260 nm. The alumina membranes were coated by atomic 
layer deposition (ALD), that consists of the sequential 
deposition of thin layers from two different vapor-phase 
reactants, into a ferromagnetic Ni layer. We used Nickel- 
tocene (NiCp2) and O3 as reactants for the deposition of 
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a thin oxide layer, which was reduced after the ALD pro- 
cess into a magnetic layer by annealing at 400 °C under 
Ar + 5% H2 atmosphere. The deposition temperature 
was between 270 °C and 330 °C with deposition rates of 
0.22 — 0.3 A/cycle. Details about the preparation method 
can be found elsewhere. [lH 

For SEM and TEM investigations, Ni nanotubes were 
deposited into the alumina membrane between two lay- 
ers of Ti02. The Ti02 layers were also obtained by ALD 
[1^ and used for adding a higher stability against oxi- 
dation to the nanotubes and for preventing their damag- 
ing in the etching process. For TEM measurements, the 
Ti02/Ni/Ti02 tubes were released by etching the mem- 
brane in 1 M NaOH and washing several times with puri- 
fied water. The magnetic properties of the Ni nanotubes 
were measured by a superconducting quantum interfer- 
ence device (SQUID). The Ni layer deposited on the top 
surface of the membrane was, in all cases, removed by 
ion milling. 
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FIG. 3: Hysteresis cycles for Ni nanotubes as a function of 
the tube diameter. The magnetic field is applied parallel to 
the tube axis. The Ni layer thickness is the same in all cases 
(11 - 12 nm). 



III. MODEL AND DISCUSSION 
A. Magnetic field applied parallel to the long axis 

As pointed out by Landeros et al. 0] , the curling rever- 
sal mode is the dominant magnetization reversal process 
in magnetic nanotubes. The rnagnetization curling mode 
was proposed by Frei et al. and has been used to 
investigate the magnetic switching of films |15| and par- 
ticles with different geometries, like spheres, prolate 
ellipsoids [3, [IB] and cylinders. However, for sim- 
plicity, expressions for the nucleation field obtained using 
infinite cylinders are used. 

In the case of a magnetic field applied parallel to the 
long axis of an infinite tube, we present an analytical 
approach to the nucleation field obtained from a Ritz 
model. Calculations are shown in the Appendix and lead 
us to write 
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FIG. 2; SEM (a) and TEM (b) images for Ti02/Ni/Ti02 
nanotubes of 180 nm diameter and a Ni layer thickness of 

11 — 12 nm. The scale bar indicates in both cases 150 nm. 

Figures 2(a) and (b) present typical SEM and TEM 
images for Ti02/Ni/Ti02 nanotubes with a diameter of 
around 180 nm and a Ni layer thickness of 11 — 12 nm. 
Figure 3 shows the hysteresis cycles for three different 
samples with diameters of 180, 220 and 260 nm; pore 
length of around 5 /xm and Ni layer thickness of 11 — 

12 nm. The measured coercivities for these dimensions 
are, in all cases, around 200 Oe (1 Oe = 10'^/47r Am~^), 
higher than the coercivity of bulk Ni (around 0.7 Oe for 
Ni). 131] 



where = y^2A/ ^iqAIq , AIq is the saturation magne- 
tization, A is the stiffness constant of the magnetic ma- 
terial, K is the magnetocrystalline anisotropy constant 
and 



(14 - 13/32 + 5/34) 
3(11 + 11/32-7/34-^/36) 
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Equation ([T]) has been previously obtained by Chang et 
al. [3] starting from the Brown's equations. They ob- 
tained a{l3) = with q satisfying 
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Here Jp (z) and Yp (z) are Bessel functions of the first 
and second kind, respectively. Equation ^ has an infi- 
nite number of solutions, out of which only the one with 
the smallest nucleation field has to be considered. 
Therefore, the nucleation field depends on a(/3), which 
is related to the internal and external radii of the tube. 
Figure 4 illustrate a as a function of /3, obtained numer- 
ically from Brown's equations, Eq. ([3]), and by means 
of Eq. ^ using our analytical approach (Ritz method). 
We observe that both results are similar, showing a per- 
fect agreement for nanowires (/? = 0) and very narrow 
nanotubes (/? — *■ 1). This behavior lead us to simply use 
the analytical expression, Eq. ([2]), to calculate the nu- 
cleation field of a magnetic nanotube. Therefore, a{(3) is 
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FIG. 4; a as a function of /3 from Eq. Q (dashed line) and 
by means of our analytical approach, Eq. ([2J (solid line). 

a decreasing function of the aspect ratio /3 — a/ R of the 
tube. It is clear that, for /3 = 0, i.e., for an infinite cylin- 
der (or nanowire), a (0) = I.OStt, as previously calculated 
by Shtrikman et al. ^M, 

As pointed by Aharoni, f2^ for a prolate spheroid with 
Of) = 0, a jump of the magnetization at, or near, the curl- 
ing nucleation field occurs. Therefore, the coercivity is 
quite close to the absolute value of the nucleation field. 
Then, we assumed here that —H^ is a good approxi- 
mation to the coercivity, : when the reversal is by 
curling, as in other works considering infinite cylinders. 
[lEIIll 

Now we investigate the validity of Eqs. (H]) and ^ 
by calculating the coercivity for different samples exper- 
imentally investigated. Table I summarizes the geomet- 
rical parameters of the arrays, measured He, and calcu- 
lated ■ In our calculations we used = 4.85 x 10^ 
A/m and if = 4.5 x 10^ J/m'"^, both taken from Ref. [22] 
at room temperature. In the same reference it is pointed 
that A ranges for any material from 1 to 2 x 10~^^ J/m. 
However, by means of field-dependent elastic small-angle 
neutron scattering (SANS), the exchange-stiffness con- 
stant A for Ni was determined by Michels et al. |i23| At 
ambient temperature A — 7.6 ± 0.3 x 10~^^ J/m was re- 
ported for nanocrystalline samples while 9.2 ±0.2 x I0~^^ 



J/m was obtained at T = 5 K. For our calculations we 
choose ^ = 7.6 X 10~^^ J/m. 



TABLE L Parameters for different Ni nanotube arrays with 
different tube diameters {2R) and 5 /im tube length. He rep- 
resents the experimentally measured coercive field. 



2R (nm) 


/3 


He (Oe) 


HY (Oe) 


180 


0.88 


200 


231 


220 


0.90 


200 


215 


260 


0.91 


200 


206 



In the measured samples, small variations of around 5 
Oe for the coercive field were detected. However, we con- 
sider these variations to be within the range of measure- 
ment error. From the work of AlMawlawi et al [2j| small 
variations of the coercivity are observed for nanowires as 
a function of the aspect ratio for aspect ratios higher than 
20. In order to observe small variations in the coercivity, 
very well geometrically characterized samples need to be 
measured with well controlled inter-element interactions. 
In our samples the center to center distance is fixed, and 
then the strength of interactions is also different from 
one sample to other, making a direct comparison difh- 
cult. The computed values for the coercivity are larger 
than the experimental data. We ascribe such differences 
between calculations and experimental results to the in- 
teraction of each tube with the stray field produce by 
the array. This field originated in the effective antifer- 
romagnetic coupling between neighboring tubes, which 
reduces the coercive field as previously demonstrated in 
nanowires. L25i,i26j As /? 1 the contribution of the stray 
field diminishes and a better agreement between theory 
and experiment is obtained. Besides, a fully realistic ap- 
proach needs to consider a finite nanotube, making much 
more complex the calculations and the expression for the 
energy. Therefore, the small discrepancy between exper- 
iments and model can be regarded as the result of our 
models simplification. 



B. Angular dependence of the coercivity 

We now proceed to investigate the angular dependence 
of the coercivity for magnetic nanotubes. We calculate 
the coercive field Hj: assuming each of the previously 
mentioned reversal mechanisms, k = C (coherent) or V 
(curling). 

Coherent-mode rotation (C) 

The angular dependence of the nucleation for a co- 
herent magnetization reversal was calculated by Stoner- 
Wohlfarth 27] and gives 
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where t = tan^ {6q) and corresponds to the demag- 
netizing factor of the eUipsoid along z. For a tube, the 
demagnetizing factor can be calculated from Escrig et al. 
[|, so that 

In the Stoner-Wolhfarth model [23| the nucleation field, 
, does not represents the coercivity, , in all cases. 
However, from the discussion on p 21 [27], the coercivity 
can be written as 

^ \ 2|i/,<r(^o = 7r/4)|-|i/^| 7r/4<0o<^/2 
Curling-mode rotation (V) 

The angular dependence of the curling nucleation field 
in a finite prolate spheroid was obtained by Aharoni. [20| 
By extending the expression for the switching field to take 
into account the internal radii of tubes, we obtain 

y (iv. - ^) ' sin^ eo + (iv. - ^) \J~e, 

(5) 

Figure 5 illustrates the coercive field as a function of 
9o for an infinite nickel nanotube. Dashed curves repre- 
sent the coercivity of a nanotube due to a reversal curl- 
ing mode. The cutoff of the curves corresponds to the 
transition angles, Oq , at which a coherent reversal mode 
appears. The most remarkable feature of these curves is 
that the general shape for nanotubes is similar to the one 
for nanowires. Differences, of course, are far from being 
negligible, and the internal radii must be taken into ac- 
count in any proper analysis of experimental data. For 
example, using Eq. ([5]) we found that for a very thin nan- 
otube {/3 — 0.9) with R = 50 nm, the coercivity is almost 
the same as in a nanowire with R — 100 nm. We also 
observe in this figure that, for Oq > 60°, a small uncer- 
tainly in the measurement of can cause large changes 
in the coercivity. 

We now investigate the dependence of the coercivity 
as a function of R. We illustrate our results with trajec- 
tories of the transition angle, 9q, in Fig. 6 for Ni. Each 
line separates the coherent reversal mode (upper) from 
the curling reversal mode (lower) . Results for nickel also 
represent iron oxide tubes because of the similar magnetic 
parameters of both materials. In the considered range of 
parameters, we observe that an increase of the external 
radii, R, or f3 (see Fig. 7) results in an increase of the 
transition angle, 9q, enhancing the region of stability of 
the curling reversal mode. However, the dependence of 
the coercivity on R is stronger than on (3. 




FIG. 5: Angular dependence of the coercivity. He, in an in- 
finite nickel nanotube with R = 50 nm. The solid line rep- 
resents the coherent rotation, the dashed line corresponds to 
the curling mode for /3 = 0, and the dotted line illustrates our 
results for the curling mode in a tube defined hy P — 0.9. 




R/nm 

FIG. 6; Trajectories of the transition angle, as a function 
of R, for 13 = (solid line) and (3 — 0.9 (dashed fine). 



IV. CONCLUSIONS 

In conclusion, by means of theoretical studies and ex- 
perimental measurements, we have investigated the co- 
ercivity in magnetic nanotubes. We have obtained Ni 
nanotubes by atomic layer deposition into alumina mem- 
branes. ALD proves to be a powerful technique, which al- 
lows us to have a very precise control of layer growth. We 
have also derived an analytical expression that allows one 
to obtain the coercivity when a magnetic field is applied 
parallel to the tube axis. Good agreement between the 
measured magnetic properties of Ni nanotubes and the- 
oretical calculations is obtained. Finally, this calculation 
has been extended to the case of an angular dependence 
of the coercivity, where a transition from curling-mode 
rotation to coherent-mode rotation has been observed. 
However, further experimental work remains to be done 
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which satisfies e^m^ + e^m^ + e^m^ = 1 + (e'') , with e 
an infinitesimal parameter. 

The total energy density E is generally given by 
the sum of four terms corresponding to the magneto- 
static -E^p, the exchange E'^, the magnetocrystalline 
anisotropy E^, and the Zeeman E^ (resulting from 
the interaction between M and an external field Ha), 
can be calculated using the well known continuum the- 
ory of ferromagnetism. [l^ Because in this case we do 
not have charges in the surface, the contribution from 
the magnetostatic energy density results equal to zero, 
^dip = 0- The exchange energy density is given by 

EY^ = (^/s)/S(Vem0^dw, with s = ttR'^ {l - f3^) 
and TUi = Mi/Mo [i ~ x, y, z). Thus, we obtain 



FIG. 7: Trajectories of the transition angle, 6q , as a function 
of /3, for _R = 50 nm (solid line) and R = 100 nm (dashed 
line) . 
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(14 - 13/3^ + 5/3^ 



in order to observe this transition. 



The magnetocrystalline anisotropy energy density can be 
written as E^ = [K/ s) J {e'^ml + s^w^) dv, so that 
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APPENDIX: NUCLEATION FIELD FOR AN 
INFINITE NANOTUBE FROM A RITZ MODEL 



We use the term curling mode here not in reference 
to an eigenfunction of Brown's equation. In our case 
we replace the spatial dependence of the curling eigen- 
function by a Ritz model that approximates the curling 
eigenmode, which turns out to be quite simple. We use 
the Ritz model previously used by Ishii et al. for infinite 
cylinders. [2l| We assume the following model for the 
magnetization 
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Finally we consider the Zeeman energy density, which is 
given by E^ — — (fioHaMo/s) J errizdv. Thus, we obtain 



1-- (11 + 11/32- 7/34 + /3«) 



Now we are able to obtain the total energy density E^ . 
The second variation in the magnetic energy density E^ 
with respect to a small deviation em from z, where z is 
the unit vector along the magnetization, must be positive 
at the equilibrium state and zero at nucleation. There- 
fore, 
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where 



(14 - 13/32 -I- 5/3^) 
3 (11 + 11/32- 7/34 + /36) 
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